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Abstract 

In the present work zero of dispersive function from Stokes' second problem 
are investigated. Stokes' second problem is a problem about behaviour of the 
rarefied gas filling half-space. A plane, limiting half-space, makes harmonious 
oscillations in the plane. The linearization kinetic ellipsoidal the statistical 
equation with parametre is used. The factorization formula of dispersion function 
is proved. By means of the factorization formula zero of dispersion function in 
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an explicit form and their research is carried out. Dependence on dimensionless 
quantity collision frequency of a plane limiting gas and on parametre equation 
are investigated. 

Keywords: Stokes' second problem, kinetic ellipsoidal statistical equation, 
separation of variables, zero of dispersion function, eigen solutions, continuous 
and discrete spectrum, factorization of dispersion function. 



Introduction 

Problem of generation of shift waves by fluctuating plate or Stokes' 
second problem for the continuous matter has been formulated in the 
middle of XIX century [T]. Then, after introduction by Maxwell and 
Boltzmann the kinetic equations, Stokes' second problem began to be 
studied for the rarefied gas. 

The detailed history of this problem is stated in our works |2]-[I]- 
In these works the Stokes' second problem was solved analytically. The 
known kinetic BGK-equation (Bhatnagar, Gross, Krook) thus was used. 

In work [5] zero of the dispersion function corresponding to BGK- 
equation have been investigated. The range of values of frequencies of 
oscillations of the plane is found out, in which to within one percent exact 
expression of zero of dispersion function can replace it by asymptotic 
approximation taken by by means of expansion in asymptotic series of 
dispersion function in neighbourhood of infinitely remote point. 

Necessity of this replacement of exact expression of zero of dispersion 
function by it asymptotic representation speaks that fact, that calculation 
of macrocharacteristics of this problem demands variety calculation of 
values of composite functions in zero of dispersion function. 

Let's underline, that this problem draws to itself wide attention of 
many authors (see, for example, [Sj-fll]). It has been solved by the 
numerical and (or) by approximate methods. In works [3, [8] Stokes' 
second problem has been successfully applied in nanotechnology. 

In work [12] Stokes' second problem has been solved analytically with 
help of ellipsoidal statistical equation. 
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Methods of calculation of zero of dispersion functions for the transport 
equations of neutrons have been put in pawn in work [13]. Then these 
methods were applied and developed for various problems in works 



In the present work in explicit form eigen solutions of ellipsoidal 
statistical equation are presented. These solutions correspond to discrete 
spectrum. For this purpose in the explicit form zero of the dispersion 
equation are found. For finding of zero it is used factorization of dispersion 
function. For this proved boundary value Riemann problem from the 
theory of functions of complex variable. Coefficient of Riemann problem 
is the relation of boundary values of dispersion function from above and 
from below on the real semiaxis. 

At small values of frequency of oscillations of a plane limiting rarefied 
gas the simple asymptotic formula for calculation of zero of dispersion 
function is found. Graphic research of modules of zero of dispersion 
function is carried out, and also the real and imaginary parts calculated 
by exact and asymptotic formulas. The function of errors representing 
a relative deviation of the module of asymptotic representations of zero 
from the module of its exact representation is entered. 

The interval of values of frequency of oscillations of a plane, in which 
the value of function of errors does not exceed one percent is found out. 

In item 2 of the present work properties of dispersion function are 
studied in complex plane. In item 3 exact formulas for calculation of zero 
of dispersion function are deduced. Properties of these zero as functions 
of dimensionless frequency of oscillations of the plane limiting rarefied 
gas also are investigated. 

In work [12] Stokes' second problem has been analytically solved with 
using ellipsoidal statistical equation, which after linearization and of 
some simplifications it is reduced to the equation 

oo 

fi— h zoh(xi, a) = — -= / exp(— // 2 )(1 — a[i[i')h(xi, /i')d/i' ', (1.1) 

oxi J 

— oo 
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where 

Z{) = 1 — 

x\ is the dimensionless coordinate, x\ = x/l, I is the mean free path of 
gaseous molecules, Q = lot = — , r = 1/V, is the frequency of collision 
of gaseous molecules, to is the oscillation frequency of plates, limitting 
half-space filling of rarefied gas, a is the number parametre of problem, 
^ a ^ 1. 



1 Statement problem 

Let rarefied one-atomic gas fills half-space x > over plane solid 
surface, laying in the plane x = 0. Surface (?/, z) makes harmonical 
oscillations lengthwise an axes y under the law u s (t) = uoe~ tujt . 

We will be linearize the distribution function of gaseous molecules 
believing 

f(x,t,v) = f (v)(l + (p(x,t,v)). 

Here fo{v) = n(/3/7r) 3 / 2 exp(— (3v 2 ) is the absolute Maxwellian, f3 = 
m/(2/cT), k is the Boltzmann constant, T is the temperature of gas, n 
is the concentration (number density) of gas, m is the mass of molecule 
of gas. 

Let further v = 1/t is the collision frequency of gaseous molecules, 
r is the time between two consecutive collisions of molecules, u y {x 1 t) 
is the mass velocity of gas, a xy {x^t) is the component of viscous stress 
tensor, 

(J xy (x J t) = m J v x v y f(x,t,v)d 3 v. 

Concentration of gas and its temperature are considered as constants in 
linearization statement of problem. 

We introduce dimensionless velocities and parametres: dimensionless 
velocity of molecules C = v^ v ((3 = m/ (2kT)), dimensionless velocity 
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of gas U y (x, t) = y/]3u y (x, t), dimensionless time t\ = vt and dimension- 
less velocity of surface U s (t) = Uoe~ lujt , dimensionless component of 
viscous stress tensor 

P 

where Uq = yfj3uQ is the dimensionless amplitude of oscillation velocity 
of half-space border. Then the linearization ellipsoidal statistical kinetic 
equation (see, for example, [21]) (short: the ES-equation) can be written 
down in the form 

^ + C x f^ + <p(x h t h C) = 

Ot\ OX\ 

= 2C y U y (x h t 1 ) - 2aC x C y P xy (x h t 1 ). (1.1) 

Here a is the parametre of equation, and at a = 1 Prandtl number 
is true (Pr = 2/3). Let's notice, that for dimensionless time U s (t\) = 
U e- m \ 

Let's underline, that the problem about gas fluctuations resolves 
in to linearization statement. Linearization of problems it is spent on 
dimensionless mass speed U y (x\,t\) provided that \U y {x : t\)\ <C 1. This 
inequality is equivalent to the inequality 

\u y {x h ti)\ <C v T , 

where vt = 1/v/^ is the heat velocity of molecules, having an order of 
velocity of a sound. 

These quantities of dimensionless mass velocity and components of 
the viscous stress tensor through function (p are expressed as follows 

U y (x h t x ) = _L /" eM-C 2 )C y ip{x h t h C)d 3 C : (1.2) 



and 



= -j^ J ex V (-C 2 )C x C y< f(x h t h C)d 3 C. (1.3) 

Considering, that plate oscillations are considered along an axis y, 
we will search function <p in the form 

<p(x h t h C) = Cye- mh h(x h C x ). (1.4) 
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By means of (1.4) it is received the following boundary problem 




1.5) 



h(0,ii) = 2U , /i > 0, z = l-iQ, (1.6) 

/i(+oo,/i) = 0. (1.7) 

It is easy to show, that equation parametre a and Prandtl number are 
connected by equality 

2 . 2(1 - Pr) 

Pr = , from which a = — . 

2 + a' Pr 

To the correct (true) Prandtl number Pr = 2/3 is answered value of 
parametre a = 1. At a = the ellipsoidal statistical equation passes in 
to BGK-equation with Prandtl number Pr = 1, i.e. at Pr = 1 a = 0. 



2 Eigen solutions of continuous spectrum 

Separation of variables in the equation (1.5) is carried out to the 
following substitution 

h v (x h fi) = exp ( - ^pj$(77,/x), (2.1) 

where rj is the parametre of separation, or spectral parametre, general 
speaking, it is complex one. 

Substituting (2.1) in the equation (1.5) it is received the characteristic 
equation 

z Q (r] - (jl)$(t}, /d) = -=77710(77) -^a/irjn^r]), (2.2) 

where 

00 

n k(ji) = J exp(-/i /2 )^ /fc $(7/,/i / )d/i / , A; = 0,1. 
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From the equation (2.2) we find, that it is possible to present it in 
the form 



where 



(7/ - /i)$(7/,/i) = 



b = 



1 



rin Q {r])(l - 6/277), 



2.3) 



Further we will accept the following normalizing condition 

00 

n o(r/) - y exp(-/i /2 )$(r/, - z . 

—00 

Then the equation (2.3) has at 77, ji G (—00, +00) the following solution 



$(77, /i) = _L P I?(L_^) + e^Afi,)^ - /2) 

77-/2 



7T 



[2.4) 



where £(3;) is the Dirac delta function, the symbol Px^ 1 means an 
integral principal value at integration a: -1 , X(z) is the dispersion function, 
introducing by equality 



\{z) = 1 - ifi + 



e r2 (l — bzr)dr 



T — Z 



This function can be transformed to the form 



X(z) = -iQ + (1 - bz 2 )X (z), 
where Xq(z) is the known function from theory of plasma, 



DC' 



X () (z) 



1 f e T rdr 



7T J T — Z 

00 



Eigen functions (2.4) are called as eigen functions of continuous spectrum 
for the spectral parametre 77 fills continuously all real axis. 

Thus, eigen solution of the equation (1.5) look like (2.1), in which 
function $(77,72) is defined by equality (2.4). 
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On the condition of our problem we search the solution which is not 
increasing far from the wall. In this connection the spectrum of the 
boundary problem we will name positive real half-axes of parametre 77. 

Let's result Sokhotsky formulas from above and from below on the 
real axis for dispersion function 




-0.5 1 1 1 

-3 3 

X 



Fig. 1. Real part of boundary values of dispersion function A + (x), frequency is 
equal to = 0.637; curve 1 corresponds to Prandtl number Pr = 2/3; curve 2 
corresponds to BGK-equation, Pr = 1. 
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Fig. 2. Real part of boundary values of dispersion function A~(x), frequency is 
equal to Q = 0.637; curve 1 corresponds to Prandtl number Pr = 2/3; curve 2 
corresponds to BGK-equation, Pr = 1. 




Fig. 3. Imaginare part of boundary values of dispersion function A + (x); curve 1 
corresponds to Prandtl number Pr = 2/3, Q — 0.1; curve 2 corresponds to 
Pr = 2/3, Q = 1; curve 3 corresponds to BGK-equation, Pr = 1, Q = 0.637. 
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Fig. 4. Imaginare part of boundary values of dispersion function A~(x), curve 1 
corresponds to Prandtl number Pr = 2/3, Q = 0.1; curve 2 corresponds to 
Pr = 2/3, Q = 1; curve 3 corresponds to BGK-equation, Pr = 1, Q = 0.637. 

Difference of boundary values from above and from below on the real 
axis to dispersion function from here it is equal 



The singular integral in these equalities is understood in sense of the 
principal value. 

3 Structure of discrete spectrum 

Let's show, that the discrete spectrum consisting of zero of the dispersion 



A+O) - A"(/i) = 2 v ^r/ie- /i (1 - b/i 2 )i 



the half-sum of boundary values is equal 




— oo 
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equations X(z) = 0, contains two zero —770 and 770, from which are 
designated through 770 that zero, at which Re rjo > 0. 

At first we will consider the case of small values Q. Let's expand 
dispersion function in asymptotic series on negative degrees variable z 
at the vicinity of infinitely remote point 

w n ^ b 1 3b . _ . 

A(z) = -^ + --^ + ^ + ..., z^oc. (3.1) 

From expansion (3.1) it is visible, that at small values Q dispersion 
function has two complex zero differing only signs. We will replace a 
number (3.1) its partial sum 

X as (z) = -iQ + ^ 



2 2z 2 ' 

Then from the equation X as (z) = we will find asymptotic representation 
of zero of the dispersion equation 



1 + 3iVta/2z Q / 1 - iQ + 3iQa/2 
(n) = */' n(2 + a/V = f 20(l-»Q + a /2) ' ° < ° < L 

From here it is visible, that at Q — > both zero of dispersion function 
have as limit one infinitely remote point rji = 00 of two order. 

Now we investigate the case of any values Q. Further is required to 
us the function 

rM - *M - -^ + (l-^ 2 )A + (r) 
1 j " A-(r) ~ -ifl + (1 - 6r 2 )A (r)' 1 j 

Let's separate at function G(r) the real and imaginary parts. Let's 
notice, that 

6 = 61 + 262, 61 



1 + n 2 ' z 1 + n 2 ' 



Af(r) = Z(r)±w(r), s(r) = V^re^ 



1 



/(r) = l-2r 2 J 
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Now equality (3.2) can be written follows 

-in + [(1 - &1T 2 ) - lb 2 T 2 ]{l{T) + 1 S {t)) 



G(t) 



or 



where 



iQ + [(1 - &it 2 ) - ib 2 r 2 ]{l{r) - is{r)Y 
p + q-i(Vt-pi + qi) 



p- q-i(Q+p 1 + qi) 



p( T ) = (1 - 6ir 2 )/(r), q(r) = & 2 r 2 s(r), 
pi(r) = (1 - 6ir 2 )s(r), gi (r) = & 2 r 2 /(r). 
Now the function G(r) it is possible to present in the form 

G(t) = G 1 (t) + iG 2 (t), 

where 

Gl{T) = WY 2{) = WY 

9i{r)=p 2 -q 2 + n 2 -p 2 1 + ql 
g 2 {T) = 2[pp l + q{n + q l )l 

9o(t) = (p ~ qf + (fi + Pi + qif. 
These functions gj{r)(j = 0, 1, 2) will be necessary in the explicit form 

9i (t) = n 2 - [s 2 (r) - l 2 (r)][(l - hr 2 ) 2 + b 2 r% 

g 2 (r) = 2s(r){Q6 2 r 2 + /(r)[(l - hr 2 ) 2 + 6 2 r 4 ]}, 
^o(r) = Q 2 +2Q[(l-6 1 r 2 ) S (r)+6 2 r 2 /(r)] + [/ 2 (r)+ S 2 (r)][(l-6 1 r 2 ) 2 +6 2 r 4 ] 
In these equalities 

(l-6 lT ) +6 2T = . 

Thus, it is definitively received 

n A -n 2 Sl (T)-s (r) 



9\{r) 



l + Q 2 
where 



8 (t) = s z (t) - P(r), Sl (r) = s (r)(l - ar 2 ) 2 - 1, 
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and 



g 2 (r) = ^{ - oflV + /(r)[l + - ar 2 ) 2 ]} : 

l + Q 2 (l-ar 2 ) 2 



2s(r 
+ 

g (T) = n 2 + [l 2 (r) + s 2 (T 
aQr 2 



+2Q 



1 



1 + D 2 
aQr 2 



SIT 



It) 



l + nv v y 1 + Q 2 

It is possible to show by means of principle of argument similarly 
that, as it is made in [2], that number of zero of dispersion function it 
is equal 

oo oo 

7V = i/ dlnG w = i/ dlnG(r) = 



7T 



argG(r) 



7T 



argG(+oo) = 2x(G), 



i.e. to the doubled index of function G(r). 

Let's enter the angle 0(r) = argG(r), which is the principal value of 
argument, fixed in zero by condition 6(0) = 0, 



6(t) = arcctg ^ 777—^; = arcctg ^ 



Im G{t) 



(3.3) 



From equation g\(r) = we find its positive root 



Ma) 



^£iM + ^(£lM) 2 + ^ sn (r,a). 



Let's enter the allocated frequency of oscillations of the plate limiting 



gas 



Q*(a) = max Q(r, 

0<r<oo 



(3.4) 



This frequency of oscillations we will name critical. 

Similarly [2] it is possible to show, that in the case, when frequency 
plate oscillations less than critical, i.e. at ^ lu < Q*(a), the index of 
function G(t) is equal to unit. It means, that number of complex zero of 
dispersion function in a plane with cut along real axis, equally to two. 
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X 



Fig. 5. Angle 9 = 6(x), Pr = 2/3, curve 1 corresponds to Q = 0.1; curve 2 



In the case when frequency of oscillations of the plate exceeds the 
critical (to > £T(a)) the index of function G(t) is equal to zero: >c{G) = 
0. It means, that dispersion function has no zero in top and bottom of 
half-planes. In this case discrete (partial) solutions the initial equation 
(1.9) has no. 

Thus, the discrete spectrum of the characteristic equation, consisting 
of zero of dispersive function, in the case ^ ft < ft* (a) there is a set 
from two points t]q and —770. At Q > ft* (a) the discrete spectrum is the 
empty set. At ^ Q < Q*(a) decreasing eigen solution of the equation 
(1.9) looks like h m (x\, fi) = e _XlZ °/ 77o $(?7o, /i), where 



is the eigen solution of characteristic equation. 

It means, that the discrete spectrum of the considered boundary 
problem consists of one point 7/0 in the case < ft < ft* (a). 

At ft — > both zero ±770 as already it was specified above, move to 
one and same infinitely remote point. It means, that in the case ft = 



corresponds to f2 = 0.5. 
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the discrete spectrum of this problem consists of one infinitely remote 
point of frequency rate two also is attached to the continuous spectrum. 
In this case discrete (partial) solutions exactly two: 

2 

hi(x h ji) = 1, h 2 (x h fi)=xi- — — \i. 

I + a 

Let's result the table of critical frequencies depending on values of 
Prandtl number and equation parametre a according to (3.4). 
Table of critical frequency values . 



Prandtl number Pr 


Parametre a 


Critical frequency fT 


1 





0.733 


0.952 


0.1 


0.717 


0.909 


0.2 


0.717 


0.870 


0.3 


0.691 


0,833 


0.4 


0.681 


0.800 


0.5 


0.672 


0.769 


0.6 


0.662 


0.741 


0.7 


0.654 


0.714 


0.8 


0.648 


0.690 


0.9 


0.642 


2/3 


1 


0.637 



4 Factorization of dispersion function 

Here we deduce the formula representing factorization of dispersion 
function in the top and bottom half-planes, and also the formula for 
factorization boundary values of dispersion function from above and 
from below on the real axis is deduced. Such factorization it is given in 
terms of function X(z). 

At the heart of the analytical solution of boundary problems the 
kinetic theory lays the solution of the homogeneous boundary value 
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Riemann problem (see [23]) with coefficient G(fi) = A + (/i)/A (//) 

85 - " > °- 

Homogeneous boundary value Riemann problem is called also (see [23] ) 
the factorization problem of coefficient G(fi). 

The Riemann problem mean that relationship A + (/i)/A~(/i) can rep- 
lace by relation X + ([i)/X~([i). Here X ± (fi) are boundary value of 
function X(z), analytical in complex plane C and having as jumps line 
the real positive half-axes. Dispersion function have as line of jumps all 
real axis. 

We consider corresponding homogeneous boundary value Riemann 
problem 

X + (fi) = G(fi)X-(fi), fi>0, (4.1) 

where coefficient of problem G(r) is defined by equality (3.2). 

The solution of Riemann problem PniviaHa (4.1) is carried out similarly 
[3] and given by integral of Cauchy type 

X(z) = ±expV(z), (4.2) 

where x = x{G) is the index of coefficient of problem, entered in item 
3, and V(z) is understood as integral Cauchy type 

o 

Here lnG(r) = In |G(r)| +i9(r) is the principal branch of logarithm, 
fixed at zero by condition lnG(O) = 0, angle 9{t) = argG(r) is the 
principal value of argument, entered by equality (3.3). The integral (4.3) 
is more convenient to consider in the form 

oo 

= 1 / jM-^tx 

7T J T — Z 


where 

9W = ^-|ln|G(r)|, 
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or in the form 



00 



C(r)dr 



T — Z 



X = 0,1, 



where 

CM = 9(0 - ^ 

Let first = 1, i.e. n £ [0, f2*(a)). We show that for dispersion 

function A(z) everywhere in the complex plane C, excepting the real 
axis R, is true the formula 



X(z) = -X(oo)(z 2 - V l)X(z)X(-z] 



(4.4) 



Here 



A(oo) = —iQ + 

in i - 



b 

2 " 
zftPr 



1 + 



2(1 -in) 



^ Pr < 1. 



Pr ' 3 

From this formula follows that for its boundary values on R are 
carried out the following relations 

X ± ( f i) = -A(oo)(/i 2 - r /0 2 )X ± ( / i)X(-/i), fjt > 0, (4.5) 

= -A(oo)(/i 2 - 7/ 2 )X(/i)X^(-/i), /x < 0. (4.6) 
For the proof of the formula (4.4) we will enter auxiliary function 

AM 



R(z) = 



(4.7) 



-\{^-r ] l)X{z)X{-zy 

This function analytic everywhere in the complex plane, except points 
of cuts R + and R_. These points z = ±770 are removable, because at 
these points A(±r/o) = 0. 

Each point of cuts R + and R_ is removable. Really, if fi > 0, then on 
the basis of equality (4.4) and (4.7) it is had 

A+M = A_(/i) 

A(oo)(^ - nDX+i^Xi-fi) A(oc)(^ _ r / 2 )X-(/i)X(-/i) ' 
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from which R + (fi) = R (//), ji > 0. If ji < 0, then on the basis of 
equality (4.4), in which fi is replaced on — we have 

X-(-m) A-W P<U ' 
It easy to see that A + (— fi) = A~(/i), A~(— /i) = A + (^i). Hence 

x+HQ A-( M ) _ 

from which 

A+(m) A-(//) 



A(oo)(/i2 - 77 2 )X(/i)X-(-//) A(oo)(^ _ 



ii < 0, 



or 



= iT(/i), /i < 0. 

To prove equalities 



and 

R-(p) = A_(M) 



-A(oo)(p 2 - % 2 )X( M )X+(-rt' 
we notice that if poin z tends to point fi (/i < 0) from upper or below 
half-plane, then functions R + {fi) or R~(fi) are calculated according to 
the previous equalities. 

Hence, it is possible to consider this function as analytical function 
everywhere in C, and in cut points, having predetermined it on the cut 
by continuity. It is necessary to notice, that function R(z) is analytical 
everywhere in C and R(oo) = 1. Under Liouville theorem this function 
is identically constant: R(z) = 1, whence the formula (4.4) is proved. 
Formulas (4.5) and (4.6) obviously follow from the formula (4.4). 
From the formula (4.4) we will find in the explicit form the formula 
for calculation of zero of dispersion function 



z* + 



X(oo)X (z)X(-z) 
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In this formula as a point z it is convenient to take the point on the 
imaginary axis: z = iVz, N = 1, 2, • • • . Then we will receive the formula 



m (n) = J-N* + 



X(Ni) 



X(oo)X(Ni)X(-Ni) 



We will calculate both parts of equality (4.4) at the point z = i. As a 
result for zero of dispersion function it is received the following formula 



rj (Q) = + ^ exp [ - V(i) - V(-i)] 



(4.8) 



Let's consider the case of zero index: x(G) = 0, i.e. Q G (fT,+oc). 
Similarly previous are proved formulas 

\{z) = -X(oo)X(z)X(-z), Imz^O. 

A ± ( /U ) = -X{oc)X ± {fi)X{-fi), fi^O. 

A ± ( /i ) = -A(oo)X(/i)X T (-/i), fi<0. 

Research of properties of zero of dispersion function is carry out on 
to the formula (4.8). On fig. 5 it is spent comparison of modules of exact 
value of zero r]o(fl) (curve 1) and asymptotic representations zero t]q s 
(curve 2). 

Let's enter function of errors, which is the function of relative deviation 
asymptotic representations of the module of zero from the module of its 
exact representation 




On fig. 7 the behaviour of function of errors is presented as function 
of dimensionless frequencies of oscillations of the plane limiting rarefied 
gas in Stokes' second problem. From fig. 7 it is visible, that in an interval 
0^ £1 ^ 0.2 quantity of function of errors not exceeds one percent. This 
fact allows in applied questions to use asymptotic representation of zero 
of dispersion function. 

5. Conclusions 
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Fig. 6. Comparison of modules of zero 770 = 9q(Q) (curve 1) and rj = r]Q S (Q) 

(curve 2), Pr = 2/3. 

In the present work zero of dispersion function from the Stokes' second 
problem are investigated. Stokes' second problem is the problem about 
behaviour of the rarefied gas filling half-space. A plane, limitting the 
half-space, makes harmonious oscillations in own plane. It is used the 
linearization kinetic ellipsoidal statistical equation. By means of the 
solution of boundary value Riemann problem the formula of factorization 
of dispersion function is proved. By means of the factorization formula 
in explicit form there are zero of dispersion function and their research of 
dependence on quantity of dimensionless frequency of the plane limiting 
gas is carried out. The interval of values of frequency of oscillations of 
the plane 0^ £1 ^ 0.2 in which the quantity of function of errors does 
not exceed one percent is found out. 



References 



21 



G 
O 



s - 



-1.5 - 



-1.0 - 



-0.5 - 



0.0 - 




-2.0 



0.1 



0.2 



Fig. 7. Relative deviations of the module of zero |?7Q S (f2)| from |?7o(fi)|, Pr = 2/3. 
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